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Q\ ' Abstract. We present a sample of mid-infrared detected sources from the European Large Area ISO Survey 

(ELAIS) regions characterised by strong mid-IR radiation with faint near-IR and optical counterparts. These 
extreme mid-to-near-IR objects (EMNOs) are defined here by a flux ratio of /„(15^im)//„ (2.2/im) > 25. This 
population is not obvious in deeper small area ISO surveys, though it produces more than 20% of the observed 
cosmic IR background radiation (CIRB) at 15/j.m above 1 mjy. Near-future large area deep mid-IR surveys with 
the Spitzer Space Telescope, however, are bound to uncover large amounts of these objects, which we argue to 
most likely be obscured AGNs, based on SED shapes and X-ray data. Very strong dusty starbursts at z > 1 may 
also have high mid-to-near-IR flux ratios, but using the MIR/NIR and FIR/MIR ratios these may be separated. 
Most of our EMNOs appear to be ULIRGs, half are also extremely red objects (ERO). A curious case of a low 
redshift, less luminous object with a very young stellar population is also found. We predict that the simple broad 
band selection method makes EMNOs a useful window into high-redshift obscured nuclear activity and its sought 
after relation to star-formation, in a similar way that EROs have been used to define samples of high-redshift 
early type galaxies. 
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1. Introduction tory of the universe at z = 1 — 3 for systematic and sta- 
tistical studies. 

g . The broadband photometric study of galaxies has turned Luminous and u l tra luminous IR-galaxies, (U)LIRGs 

out to be an efficient method of discovering new classes ^ > 1Q n _ 1Q 12 L ) are a key popu i a tion producing 

of galaxies. For example, the extension of optical pho- mQst of thfi gteUar energy Qutput m the Universe since the 

tometry into the near-infrared region revealed a popula- recombination era (Elbaz & C esarskyG>003). On the other 

tion of sources that was notrepresented m optical surveys hand; geyeral shortcomings of the generally successful and 

(Elston, Rieke & RiekeEH^. These objects were called appealing hierarchical scenarios of galaxy formation may 

"Extremely Red Objects" (EROs) because of their very be due tQ the poQr unde r S t a nding of the AGN/galaxy 

red optical-infrared colours (R - K > 5, / - K > 4), and formation and starburst interplay (Granato et al. [2004). 

have become a useful tool in constraining galaxy forma- Indeed; extreme IR _ galaxies a i ways see m to have AGN 

tion and evolution models in the redshift range = 1-3. contributions in varying deg rees (e.g. Farrah et al. M 

The Infrared Space Observatory {ISO) opened up Serjeant et al. I2003f> . 
wavelengths for large scale gala xy sur veys further in the In tMg p aper we prese nt a sample of faint ISO de- 
infrared (e.g. Franceschini et al.|200I|. However, only the tected mid _ IR ga l ax i es which we argue to belong to these 
shorthly available wealth of mid- and far-infrared photo- crucial classes of ob j ects needed to be understood before 
metric data obtained with the Spitzer Space Telescope ga l ax y evolution and structure formation in the Universe 
(Werner et al. |2004| will truly uncover the obscured his- can be fully modelled. We show that a selection made 
by high mid-IR to near-IR flux ratio will uncover specif- 
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ponent growing with the afore mentioned ratio. Nearly 
all the sources found in this study are above the ULIRG 
limit. We discuss sub-classes of these bright mid-IR ob- 
jects to aid the classification and study of IR-galaxies un- 
doubtedly soon to be cataloged by Spitzer in plentitude. 
Cosmology used throughout is: Vt m = 0.3, fl\ = 0.7 and 
H = 70 kms^Mpc" 1 . 

2. The EMNO criterion 

We search for objects with a flux ratio 
/ l/ (15/xm)// 1/ (2.2^m) = fis/fg; > 25, and dub these 
extreme mid-to-near- IR objects, EMNOs. The limit was 
chosen because the very strongest starburst galaxies, such 
as Arp220, are expected to have approximately this flux 
ratio (e.g. Charmandaris et al. 2002), while AGN related 
phenomena produce steeper NIR-MIR SEDs (Haas et al. 
127)1731 Prouton et al. 127)01)) . 

The adopted criterion is tested in Fig. fusing ISO data 
and models. In the upper panels we plot the MIR/NIR 
flux ratio against the r — K colour. The small crosses in 
the left panel are (extragalactic) sources from European 
Large Area ISO Survey (ELAIS; Oliver et al. l27)07)j> 
taken from the multiwavelength ELAIS Band-Merged 
Catalogue 1 (Rowan- Robinson et al. l2004f) . They are shown 
to mostly consist of "normal" starforming galaxies (see 
also Pozzi et al. 2003 Vaisanen et al. I2002f) and fall into 
the lower left quadrant. The flux limit of ELAIS at 15^m 
is approximately 0.7 mJy, and K w 18 and r *=s 24 are the 
limits in the other bands. 

We also checked deeper ISO surveys, with sensitivity 
limits varying between ks 0.1 — 0.4 mJy, in the HDF- 
N, HDF-S, Lockman Hole, and CFRS fields (data from 
Aussel et al. IT533I Mann et al. 127)1721 Fadda et al. 127)011 
Flores et al. H999f) . Data with appropriate bands available 
are overplotted in the upper right panel of Fig.^ Though 
these deeper and small area surveys on average pick up 
slightly redder sources than ELAIS, they do also clearly 
concentrate in the same lower left region. The only popu- 
lation with /15//K > 25 were found from the Fadda et al. 
I|2002l see also Franceschini et al. 2002) X-ray selected 
Lockman data-set of ISOCAM counterparts to XMM- 
Newton sources. All of the Fadda et al. f 15/ fit > 25 
sources are classified as obscured AGN (typically with 
Nh > 10 22 cm~ 2 ) using X-ray characteristics, though 6 
out of 7 of them remained unclassified in the optical. 

Galaxy model colours, plotted in the lower panels of 
Fig.[JJ are calculated using SEDs from the GRASIL model 
(Silva et al. [1998). As seen, normal starforming galaxies 
always have fis/fx < 10- IR-galaxies with very strong 
PAH emission such as Arp220, have larger flux ratios, 
upto our EMNO criterion. They are rare in the local 
Universe. ULIRGs at higher rcdshift are expected to be 
picked up more readily. However, K-correction (combined 
with strong silicate absorption) hampers their detection 

1 http:/ /astro. imperial. ac.uk/Elais/Data_release/ 



at z ~ 0.5 resulting in an apparent gap in the MIR/NIR 
distribution in a flux-limited sample. 

On the other hand, in sources where AGN activity 
heats the dust, MIR/NIR ratio is typically large at any 
redshift. It is important to stress that the EMNO selection 
picks out specifically obscured AGN. Obscured nuclei with 
Njj values in excess of 10 22 cm -2 , are always EMNOs, 
often with /ib//k > 100 and the ratio increases with in- 
creasing Njj (calculated from SEDs in Silva et al. 12.004), 
while e.g. Syl nuclei are never EMNOs (Fig. ^ dash-dot 
curve). We note that in the parameterization of SEDs of 
Seyfert nuclei between 1 and 16/xm by Alonso-Herrero et 
al. (2003), our EMNO limit corresponds to their spectral 
index am > 1.7, where their Sy2s are found. 

Inclusion of the host galaxy properties complicates the 
models. In Fig. the higher plotted 'AGN-2' model is a 
'custom scaled' composite of Arp220 SED and the nuclear 
SED of NGC1068 (latter taken from Silva et al. l27M|l . 
resulting in an EMNO. The lower solid curve is calcu- 
lated from composite type-2 QSO + host galaxy models 
of Silva et al. H2004[) . using a combination giving a max- 
imal /15/ fjc: AGN X-ray luminosity Lx > 10 44 ergs s _1 
and 10 23 < N H < 10 24 . This implies that most local 
AGN-hosting galaxies, where the SED templates were de- 
rived from, are not quite extreme enough to be EMNOs. 
The luminosity of the host galaxy vs. the AGN plays an 
important role - and thus can potentially be studied - 
with EMNO selection: e.g. in the models of Treister et al. 
(2004), where an evolved elliptical is the host of all ob- 
scured AGN, there are differences by factors of several in 
/15/ Jk compared to models used above with equivalent 
choices of Lx and Njj. However, the general characteris- 
tics of the models are similar, i.e. with host+AGN com- 
bination models EMNOs are expected only from sources 
with strong AGN contribution, Lx > 10 43 ergs s _1 , with 
the highest fis/fa ratios at z > 1 produced by sources 
with Nh close to values of 10 23 cm~ 2 . 

When using the Spitzer MIPS/IRAC bands, Arp220 
-like ultraluminous starbursts have /_4//3.4 < 35 (ex- 
cept locally), and our EMNO definition would then be 
j 24/ '/3.4 > 35, corresponding to Spitzer colour of [3.6 — 
24] > 7.8. Indeed, a handful of such sources are evident 
in early Spitzer results (e.g. Lonsdale et al. 2004, Chary 
et al. 121701 Yan et al. 121101 Le Floc'h et al. E_U04) - they 
are tentatively interpreted as being luminous IR-galaxies 
where an AGN cannot be ruled out with the available 
data. The mentioned MIR/NIR flux ratio does not signif- 
icantly change with the used IRAC band. With if-band, 
/24//.R" > 60 would be appropriate. It is noteworthy that 
when using the longer wavelength IRAC bands the ex- 
pected AGN-2 flux ratio becomes increasingly similar to 
strong starbursts, and thus if -band or 3.4^tm is best if the 
goal is to separate AGN from starformation. 

3. Observational data and the sample 

Our sample of near- and mid-IR detected galaxies pre- 
sented here are from the ELAIS fields Nl and N2. 
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Fig. 1. Left panel: Upper panels: The r — K colour vs. the 15/im to -ftf-band flux ratio. Solid symbols show our EMNOs, 
circles identifying IRTF EMNOs and squares other ELAIS EMNOs. Non-EMNO ('nE') IRTF and ELAIS sources are 
plotted as indicated, as well as data from HDF-S (Mann et al. 2002), CFRS (Flores et al. I1999[) . and X-ray selected 
MIR-sources from HDF-N and Lockman hole (Fadda et al. E20Q2) ). Lower panels: model colours for fis/ fx (ISO) and 
fill Ha (Spitzer), respectively, vs. redshift. Models are calculated using SEDs from GRASIL (Silva et al. I1998[) and 
Silva et al. (see text). 



The 15/zm data, taken in the ISOCAM LW3 band and 
reaching approximately 0.7 mJy, are from ELAIS 15/im 
Final analysis catalogue (Vaccari et al. 12004(1 . available 
also within the multiwavelength ELAIS Band-Merged 
Catalogue (Rowan- Robinson et al. 2004). The NIR obser- 
vations, however, are deeper than those available in the 
Band-Merged Catalogue and are described in Vaisanen 
& Johansson (2004b) in detail - essentially, we carried 
out a small-area survey to K < 20 mag targeting faint 
ISOCAM sources in approximately the range 1 to 3 mJy, 
without obvious counterparts in DSS images. A total of 
12 ISO sources were observed in this survey, and all were 
found to have K-band counterparts. Observations were 
made with the 3.0-m NASA IRTF telescope on Mauna 



Kea under photometric conditions (seeing varied between 
0.7-1.0") using the NFSCam 256 x 256 format InSb detec- 
tor array with a 0.3" pixel scale. Optical photometry in 
SDSS-like filters is taken from the Isaac Newton Telescope 
Wide Field Survey (WFS) 2 data products (McMahon et 
al-EHHl Gonzales-Solares et al. TMty . 

The MIR data is matched with the NIR-optical data- 
sets using a 4" search radius around the ISOCAM detec- 
tions. We do not expect spurious matches given that the 
probability of finding a chance K < 19 mag object within 
the search area is P w 0.03, and the ISOCAM detection 



2 http://www.ast.cam.ac.uk/ ~wfcsur/index.php 
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error circle radius (la) is better than 2" (Vaccari et al. 

12003- 

We find 6 EMNOs from the fRTF data, from the to- 
tal of 12 matches between ISO and IRTF data. One of 
them is the fainter NIR counterpart of an ISO object with 
two possible counterparts. In addition, we searched for 
EMNOs in the full published ELAIS catalogue (the ISO 
data characteristics are the same as in our EMNO survey), 
which covers a much wider area but utilises a significantly 
shallower NIR coverage: three additional sources are found 
(EMNOs 1, 2, and 6). The photometry of all the 9 sources 
is presented in Table ^ an d Fig. [21 shows the 6 EMNOs 
from our own i^-band images with 15/zm contours over- 
laid. 

4. Discussion 

4.1. Colours and colour-colour diagrams 

The MIR/NIR flux ratio vs. optical to NIR colour of our 
EMNOs are shown as solid symbols in the upper panels of 
Fig-El It is obvious that they are well separated from the 
bulk of MIR galaxies. Our non-EMNO IRTF/ISOCAM 
sources are plotted with open circles. Six of the EMNOs 
are also EROs, or very close to being EROs (defined as 
r — K > 5.5, see Vaisanen & Johansson 2004a). As is well 
known, EROs are thought to be either passively evolving 
ellipticals or highly reddened dusty starbursts. For com- 
parison, we overplot (open squares) three spectroscopi- 
cally confirmed ISO-detected EROs (Pierre et al. IMiTl 
Smith et al. 120011 Elbaz et al. 12002 b|) . None of them are 
EMNOs, although the overall SED of especially HR10 (the 
reddest r — K of the plotted three) is very similar to e.g. 
our EMNO-5 (see below). 

We note that all our EMNOs are part of the high 
fis/ fr' > 10 2 ratio population found by Gonzalez-Solares 
et al. I|2004fl and La Franca et al. Ij2004|l . mostly consisting 
of 15 /zm sources not identified in optical images. Those 
EMNOs which are also EROs have much more extreme 
MIR/optical ratios, /i 5 // r < > 10 3 . 

4.2. SEDs and redshifts 

Figure shows the observed SEDs of our EMNOs. It is 
even more obvious than from the r — K colour alone that 
there are two kinds of EMNOs, those with optically red 
colours (EMNOs 1-6) and those that are bluer (7-9). We 
discuss them separately. 

We fitted representative IR-bright GRASIL models 
(Silva et al. H^gjl . namely M82, Arp220, and the ex- 
tremely obscured star burst ERO HR10, to the optical- 
FIR SEDs of red EMNOs. Later ISO data has shown the 
HR10 IR-SED to closely resemble that of the starburst 
ULIRG Arp220 (Elbaz et al. I2002bjl but we chose here to 
retain the older GRASIL-SED to represent the intermedi- 
ate IR-luminositics between M82 and Arp220. 

We had spectroscopic redshifts for only two objects 
in our EMNO sample and for the other sources we had 



to resort to photometric redshift methods. For the three 
EMNOs with optical detections we calculated photomet- 
ric redshifts using the HYPERZ software (Bolzonella et 
al. l2TTUU|l . For EMNO-3 we found a best fit of z = 1.2±g;2 
(all errors are la) and for EMNO-7 we found a photo- 
metric redshift of z — 2.5^q'1- We also searched for sec- 
ondary redshift solutions and for EMNO-3 and EMNO-7 
we found secondary solutions that were within the la error 
bars of the primary solutions. For EMNO-5 the photomet- 
ric redshift was better constrained with a very well fitted 
Arp220-like SED at a redshift of z = O.^o l all the way to 
radio wavelengths (this is one of only two of our EMNOs 
also having a radio detection) making this source a good 
candidate for a high-redshift Arp220 ULIRG clone. For 
EMNO-5 the secondary redshift solution z sec = 1.4 had a 
confidence level that was lower by a factor of ten than the 
confidence level of the primary solution. 

For the four objects with no optical detections we se- 
lected the most 'conservative' redshift and SED combina- 
tion (i.e. starting from low-redshift M82 SEDs) fitting the 
optical and FIR upper limits when normalized in the K- 
band. Three HR10 SEDs and one M82-type SED were fit- 
ted, i.e. all highly reddened starbursts at z ~ 1.5 (Table^). 
With our current data we can conservatively estimate our 
photometric redshifts to be accurate within z — ±0.5 with 
a higher confidence given to the objects with HYPERZ fits 
than the ones with only SED fits. 

Interestingly, however, the FIR (and radio) upper lim- 
its of most red EMNOs rule out an Arp220 SED, which 
otherwise would explain most or all of the observed MIR 
radiation. This is very similar to the recent Haas et al. 
(2004) findings where sources selected by high /6.7//K 
ratio have low feo/fMiR, indicative of AGN. Our SEDs 
are clearly suggestive of hotter dust as well, as seen 
from the overplotted appropriately redshifted and arbi- 
trarily normalized nuclear SED of the well-known Sy- 
2 NGC1068 (Fig. thick line). EMNO-6 has addition- 
ally been overplotted with an example of Treister et al. 
i|2TiU4) obscured AGN model with N H = 10 23 cnT 2 and 
L x = 10 43 - 5 erg s -1 , which fits the red EMNOs, with the 
exception of EMNO-5, remarkably well. In general, the 
higher the MIR/NIR ratio is, the more likely and domi- 
nating the AGN component will be, as already suggested 
above in Section|21 For example, in the sample of Prouton 
et al. (|2004|) . the only ULIRG that would be classified as 
an EMNO is their only one where close to ~ 50% of the 
IR luminosity is provided by the AGN component. 

We do note that extreme starformation alone in ob- 
scured very young starburst regions can also produce steep 
MIR SEDs without yet contradicting FIR upper limits 
(Takeuchi et al.EOOH- At z > 1 the red EMNOs would 
truly need to be massive primeval galaxies undergoing 
their very first starburst, and we find the AGN contri- 
bution option to be the more conservative one. 

The blue EMNOs have optical colours of a less ex- 
tincted younger stellar population. The overplotted SEDs 
are those of starforming GRASIL spirals (Table EJ, which 
appear to need large contributions to MIR flux from ob- 
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Table 1. Photometry of the EMNOs. Columns (2)-(5) give optical magnitudes (Vega) in 2.4" diameter apertures. Our 
if-band (Vega) data (6) is determined in matched apertures, except sources 1, 2, and 6 which use total Jf-band mags 
from the ELAIS Band Merged Catalogue (Rowan-Robinson et al. 12004(1 . Typical photometric errors are ~ 0.10 mag 
for both optical and NIR data, and typical seeing was 1" for all data and no seeing corrections were applied to the 
magnitudes. (7)-(8) 15/zm flux and the corresponding signal-to-noise ratio. The /15/ fx flux-ratio (9) is calculated from 
total fluxes. (10) is the redshift: three decimals signal a spectroscopic redshift, one decimal a photometric determination 
(optical-NIR data) with HYPERZ la errors, and a value in parenthesis is a best-guess based on the over-all SED with 
typical errors of z — ±0.5. Columns (11)-(12) tabulate the SFR and total IR luminosity assuming the (10) redshift 
and SED given in (13). The optical-NIR counterpart of EMNO-3 is one out of two possible sources (see Fig. the 
brighter one is a non-EMNO QSO at z — 1.15. 
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Fig. 2. The 15 /Ltm contours overlaid on IRTF If -band images. Contours designate 2,3,5,7,10,15, and 20a- SNR levels. 



scured central engines. Accordingly, EMNO-9 is a spec- 
troscopically confirmed (z = 2.610) moderately reddened 
QSO with very large X-ray absorbtion (Njj ~ 3 x 
10 23 cm" 2 ; Willott et al. EUOlfl object N2_25 therein). All 
of our other EMNOs unfortunately fall out of the X-ray 
coverage of ELAIS, but the optical to MIR characteristics 
of EMNO-7, with a photometric redshift of z = 2-5±o% 
(HYPERZ la errors, X-band not used since there is likely 
AGN contribution), look very much like that of EMNO-9. 

EMNO-8 is a peculiar case, with a spectroscopic red- 
shift of only z = 0.1882, making its luminosity lower than 
the other EMNOs by two orders of magnitude. Optical 
and NIR sizes indicate a dwarf, r^i ~ 2 kpc. Interestingly, 
its SED follows well that of SBS 0335-052 (Dale et al. 
2001J an extremely metal poor dwarf galaxy, whose sur- 
prisingly strong mid-IR flux is interpreted to result from 
hot dust heated in an extremely intense interstellar radi- 
ation field around a few Myr old starforming region. Our 
EMNO-8 is approximately 10 times more luminous, and 
whether its properties are a scaled up version of SBS 0335- 
052 or whether there is some AGN activity involved, re- 



mains to be studied. Puzzlingly, the ELAIS Band-Merged 
Catalogue (Rowan-Robinson et al. 12004(1 does not men- 
tion any signs of starburst or AGN signatures in its opti- 
cal spectrum. Finally, we note that an r — K > 6.5 ERO is 
just within the 3cr error circle of the ISOCAM position - if 
it were the true counterpart, the SED and interpretation 
would be identical to the red EMNOs. 



4.3. Star-formation rates 

Due to their extreme nature, star-formation rates (SFR) 
based on indicators used in the literature cannot be con- 
fidently applied to our EMNOs. The following should 
therefore be regarded as an attempt to obtain indica- 
tive levels of SFR for the EMNOs. First, total IR lu- 
minosity is calculated from the best fit GRASIL SEDs 
of each object and by scaling the SED to observed K- 
band luminosity. Following Mann et al. 1(2002(1 we then de- 
rive SFR [M Q /yr] = L m (3 - 1000^m)/[2.3 x 10 36 W] 105 . 
These are the Lir and SFR values listed in Table ^ 
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Fig. 3. Optical-to-FIR SEDs of our EMNOs. Solid circles show detections in g', r', i', Z, K and 6.7, 15, 90, 175/z bands, 
and the open circles the corresponding upper limits. GRASIL models (thin line, Silva ct al. 1998) given in Tableware 
overplotted, as well as the nuclear SED of NGC1068 (thick line, Silva et al.P004). EMNO-1 is overplotted also with a 
3Gyr elliptical at z — 2.5 (dashed), and EMNO-6 with an obscured AGN model from Treister et al. I|2004l see text). 
The observed SED of SBS0335-052 is shown as squares with EMNO-8. For reference, the lower right panel shows a 20 
Myr SED from Bruzual & Chariot iflM^ and a 1000 K blackbody. 



SFRs using SEDs of M82 and Arp220 bracket values of 
- 200- 1000 M /yr. 

Similarly, normal late type spiral SEDs approximately 
fitting the blue EMNOs result in SFRs of the order of 
~ 100MQ/yr for the high redshift ones, and < 1 M /yr 
for EMNO-8. 

Had we scaled the model SEDs to the detected 
ISOCAM fluxes, the SFRs would have been higher by fac- 
tors of 2-5. Similarly, if we directly used the 15/im flux 
as a SFR indicator as calibrated by Elbaz et al. ( 2002a 
see their Fig. 5), our red EMNOs would have very large 
SFR values at a level of 5000 Mq / yr, and the high- z blue 
EMNOs 3-5 times higher than these. This estimate is of 
course incorrect if, as suggested by the SED shapes, most 
of the 15 /xm flux density comes from AGN related hot 
dust emission rather than from star formation. However, 
the more conservative values obtained above might not be 
any more accurate, EMNOs being very extreme objects 
with regard to their SED shape in the near to mid-IR. 

Finally, it is important to note that also fully evolved 
ellipticals at z ~l-2 satisfy very well our optical-NIR 
colour limits. In this case all the mid-IR flux would have 
to be obscured AGN originated - to illustrate this, the 
SEDs of EMNOs 1 and 6 are overplotted (dashed lines) 
with a GRASIL elliptical, and AGN+elliptical model from 
Treister et al. 1|20U4|I . respectively. Nearly all red EMNOs 
are consistent with this old elliptical host + AGN possi- 
bility, and thus the SFR lower limit for 5 out of 6 of these 
is essentially zero. From the elliptical GRASIL model we 
estimate that if the if-band light of the red EMNOs is 



assumed to be totally of stellar origin, their stellar masses 
are of the order of 0.5-1T0 12 Mq, implying black hole 
masses of > 10 9 M (Marconi & Hunt l2UU3|) . 

4.4. Number counts and contribution to the 
extragalactic background 

Proper number densities cannot be derived because of a 
targeted and discontinuous survey area. However, we es- 
timate a definite lower limit on the number of mid-IR 
selected EMNOs in the following way. In our IRTF follow- 
up 12 ISO sources (with wl-3 mJy) were observed, 6 of 
which were EMNOs (and 3 of which red EMNOs). The 
ISO sources were selected because of their faint or non- 
existent counterparts in DSS, thus obviously biasing us 
to large MIR to optical-NIR flux ratios. These ISOCAM 
sources were spread around an area of 0.15 deg 2 , which 
altogether contains 30 15/xm sources in the range 1-3 mJy, 
including the ones observed. Thus, assuming that none 
of the non-observed MIR sources are EMNOs, EMNOs 
make up 20% of all 15/im sources at ~ 1 mJy. If the tar- 
geted 12 sources are representative, on the other hand, the 
EMNO fraction could be even 50%. In contrast, less than 
1% of the general ELAIS catalogue sources (all fluxes) are 
EMNOs (Fig.[IJ. In the selected deeper (in MIR) ISO sur- 
veys (Fig.[TJ, < 10% of the sources are EMNOs - however, 
these surveys have widely varying depths of NIR coverage. 

In terms of total flux density, all the ISOCAM sources 
brighter than 1 mJy (excluding a bright star) in the IRTF 
survey area amount to rj 0.3 nW m~ 2 sr _1 , i.e. 10% of the 



Johansson, Vaisanen & Vaccari: EMNOs: obscured AGN and dusty starbursts 



7 



total CIRB at 15 ^m, according to Elbaz et al. (2002a), 
totally consistent with values therein. EMNOs contribute 
22% of the flux of all the ISOCAM sources in the area, 
which thus is the minimum EMNO contribution to CIRB 
at this flux level. As above, if the 6 EMNOs on the other 
hand were representative of the area, up to half of the 
CIRB is contributed by EMNOs. Fadda et al. lj!2UU2l) find a 
17±7% contribution by AGN to MIR background at these 
flux levels. As we have argued, EMNOs are likely to be 
dominated by AGN activity. Assuming full EMNO/ AGN 
correspondence, our value for AGN contribution is thus 
consistent with theirs. If the fraction of CIRB contribution 
using our method turns out to be larger, up to 50% as is 
possible, this could suggest that X-ray surveys are indeed 
missing the most obscured AGN, and MIR surveys could 
help in this regard (see e.g. Haas et al. 120041 Gandhi & 
Fabian 2003); Compton thick AGN would be obscured 
even in hard X-rays while their nuclear SED in NIR-MIR 
should still have EMNO characteristics. At the moment 
we cannot check this speculation because of lacking X-ray 
observations of the sources. 

5. Conclusions 

We have presented a population of sources characterized 
by very strong MIR flux compared to their optical and 
NIR properties. The sources have been dubbed EMNOs, 
and are argued to host an obscured AGN based on: i) 
Their spectral shape in the near to far-IR region which 
shows a high MIR/NIR flux ratio and in most cases a 
low FIR/MIR ratio, ii) The high-redshift spectroscopically 
confirmed EMNO is shown to be a heavily obscured quasar 
by X-ray data (others do not have X-ray coverage), iii) All 
the objects found from several other ISO surveys which 
satisfy our EMNO criterion are very hard X-ray sources 
suggesting obscured AGN. 

Those of our EMNOs which are red in their optical- 
NIR colours are well understood as z = 1 — 1.5 (ul- 
tra)luminous dusty galaxies based on their SEDs and pho- 
tometric redshifts. Their mid-IR flux density is not ex- 
plained, however, without a dominant fraction of AGN 
activity. Moreover, at least for 4 out of 6 of these cases 
are also consistent with passive early type galaxy SEDs, 
in which case all the mid and far-IR radiation would be 
AGN produced. 

The blue EMNOs appear to be higher redshift ob- 
scured quasars, according to the one spectroscopically 
confirmed case. However, their optical to FIR broadband 
colours appear to be very similar to low redshift extreme 
starbursting dwarfs, based on the similarity of our EMNO- 
8 and the well-studied metal-poor dwarf SBS 0335-052. 
Spectra, or x-ray and radio data are needed to study this 
degeneracy further. 

EMNOs contribute 20-50% of the extragalactic source 
counts and cosmic IR background radiation at the 1 mJy 
level. If further observations show them to be dominated 
by AGN, as we suspect, it could mean that a significant 
fraction of AGN activity is missed by X-ray surveys. 



We predict that the EMNO criterion will be very use- 
ful in selecting obscured AGN and sources with extremely 
strong starformation, but also in disentangling the AGN 
from starbursts since the MIR/NIR flux ratio is larger 
for the AGN with appropriately selected bands. Spitzer 
surveys using IRAC and MIPS bands (such as SWIRE; 
Lonsdale et al. I2003fl are expected to detect large num- 
bers of EMNOs. The inclusion of MIPS 70/^m data will 
be particularly useful in constraining the AGN/starburst 
fraction in EMNOs and generally in better classifying this 
new population. 
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